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Abstract 

We show that at short time aging is realized in a simple model of binary glasses. We 
find that below the transition point the off-equilibrium correlation functions and response 
functions are compatible with the relations that were originally derived Cugliandolo 
Kurchan for generalized spin glasses. 

1 Introduction 

The behaviour of an Hamiltonian system (with dissipative dynamics) approaching equilibrium 
is well understood in a mean field approach for infinite range disordered systems [|I], 0, ^. In 
the low temperature phase the correlation and response functions satisfy some simple relations 
[|T| . We present indications that they are satisfied in the case of binary glasses. 

Generally speaking in a non equilibrium system it is natural to investigate the properties 
of the correlation functions and of the response functions. Let us concentrate our attention 
on a quantity A(t), which depends on the dynamical variables x{t). Let us suppose that the 
system starts at time t = from a given initial condition and subsequently it follow the laws 
of the evolution at a given temperature T. If the initial configuration is not at equilibrium 
at the temperature T, the system will display an off-equilibrium behaviour. In many case the 
initial configuration is at equilibrium at a temperature T' > T. In this note we will consider 
only the case T' = oo. 

We can define a correlation function 



and a response function 



C{t,t^) =< A{t^)A{t + t^) > 
_ 5A{t + 



G(t, 



e=0 



where we are considering the evolution in presence of a time dependent Hamiltonian in which 
we have added the term 

' dte{t)A{t). (3) 
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The off-equilibrium fluctuation dissipation relation |I| states some properties of the corre- 
lation functions in the limit t^j going to infinity. The usual equilibrium fluctuation dissipation 
(FDT) relation tell us that 

GW ^ -/J^, (4) 

where 

G{t)= lim G{t,t^), C{t)= lim C{t,t^). (5) 
It is convenient to define the integrated response: 

R{t,t^)= f' dTG{T,t^), R{t)= lim R{t,t^), (6) 

Jo t^^oo 

i.e. the response to a field acting for a time t. 
We can also define the quantities 

C(oo) = lim R{oo) = lim R{t) = (3 (C(0) - C(oo)) . (7) 

The last equation can be naively written also as 

Rioo) = p[< > - < A>^) , (8) 

Here the brakets denote the usual equilibrium expectation value. 

In the study of off-equilibrium spin glasses systems Cugliandolo and Kurchan proposed 
that the response function and the correlation function satisfy the following relations: 

G{t, t^) ^ -PX {G{t, t^)) (9) 
The previous relation can be also written in the following form 

R{t,t^) ^ (3 / X{C)dC. (10) 

JC{t,tn,) 

The function X{C) is system dependent and its form tell us many interesting information. 

If C(oo) 7^ 0, we must distinguish two regions: (a) A short time region where X = 1 (the 
so called FDT region) and C > C(oo). (b) A large time region (usually t = 0(tu,) where 
G <G{oo) and X < 1 (the aging region) 0, ||]. 

In the simplest non trivial case, i.e. one step replica symmetry breaking, 0, ^ ^ the 
function X{G) is piecewise constant, i.e. 

X{G) =m for C < C(oo), X{G) = 1 for C > C(oo). (11) 

In all known cases in which one step replica symmetry holds, the quantity m vanishes 
linear with the temperature at small temperature. It often happens that m = 1 at T = T^. 

The previous considerations are quite general. However the function X{G) is system 
dependent and its form tell us many interesting information. Systems in which the replica 
symmetry is not broken are characterized by having m = in the formula eq. (^Ij) . 

Sometime simple aging is also assumed 0, i.e. the following the scaling relation holds 
outside the FDT region: 

Git,t^) = Gs(Pi. (12) 



Simple aging may be correct, but it is not a necessary consequence of the previous relations. 
We will see that simple aging is satisfied also in this system [0]). 
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2 The model 



The model we consider is the following. We have taken a mixture of soft particles of different 
sizes. Half of the particles are of type A, half of type B and the interaction among the particle 
is given by the Hamiltonian: 

where the radius (a) depends on the type of particles. This model has been carefully studied 
in the past [l^, |ll|]. It is known that a choice of the radius such that ctb/o"^ = 1.2 strongly 
inhibits crystallisation and the systems goes into a glassy phase when it is cooled. Using the 
same conventions of the previous investigators we consider particles of average diameter 1, 
more precisely we set 

a\ + 2{aA + asf + 0% ^ ^ ^^4^ 

Due to the simple scaling behaviour of the potential, the thermodynamic quantities depend 
only on the quantity T'^/p, T and p being respectively the temperature and the density. For 
definiteness we have taken p = 1. The model as been widely studied especially for this choice 
of the parameters. It is usual to introduce the quantity F = The glass transition is known 



to happen around F = 1.45 [10 



Our simulation are done using a Monte Carlo algorithm, which is more easy to deal with 
than molecular dynamics, if we change the temperature in an abrupt way. Each particle is 
shifted by a random amount at each step, and the size of the shift is fixed by the condition 
that the average acceptance rate of the proposal change is about .4. Particles are placed in 
a cubic box with periodic boundary conditions. We start by placing the particles at random 
and we quench the system by putting it at final temperature (i.e. infinite cooling rate). 

The main quantity on which we will concentrate our attention is the asymmetry in the 
energy (or stress): 

^ = E ^^f'^t^^^u'' (2(x. - x,f - - Vkf - {z, - z,f) = 2Ti,i - T2,2 - T3,3. (15) 
i<k ^/^l 

In other words A is a combination of the diagonal components of the stress energy tensor. If 
the particles are in a cubic symmetric box, we have that 

<A>=0. (16) 

If the box does not have a cubic symmetry the effect of the boundary disappears in the infinite 
volume limit (at fixed shape of the boundary) and we have that the stress density a vanishes 
in this limit: 

lim ^ ^ ^ = a = 0. (17) 

N 

3 Numerical simulations 

We have done simulations for various values of ranging from A^ = 18toA^ = 130. We 
report here only the results for = 130 where we have done the average over 250 samples. 
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Figure 1: On a linear scale we see the response function 5* at F = 1.8 for = 130 at different 
values of (128, 512, 2048, 8912 and 32768) as function of s = t/t^. 

those for smaller systems are qualitatively similar, and can be found in ref ||^. The evolution 
was done using the Monte Carlo method, with an acceptance rate fixed around .4. 

The limit t ^ oo may be tricky. All the previous theoretical discussion were done for an 
infinite system. In practice we need that that is greater that the time to a given power, the 
exponent of the power being system dependent. Therefore we cannot strictly take the limit 
t — * oo for finite A^. We have to study the behaviour of the system in region of time whose size 
increase with A^. Finite size corrections sometimes become much more severe at large times 
i- 

We will study the correlation C{t, tw) and the response R(t, t^). We introduce the variable 
s = t/tw According to simple aging the correlation functions should become a function of 
only s in the limit of large times. Of course the FDT region, which is located at finite t also 
when goes to infinity, is squeezed at s = 0, so that we expect that the function R becomes 
discontinuous at s = 0. Moreover the limit s — > give us information on the value of R{oo) 

lim R(t,t J = R{oo), (18) 

s— >0 

where it is understood that the limit is done always remaining in the region where t >> 1. 

We have followed a standard procedure [1^ to measure the off-equilibrium response 
function in simulations: we have kept the system in presence of an external field e up to time 
ty^ and we have removed the field just at this time. 

If e is sufficient small, we have that the stress as function of time, is related to the integrated 
response R by the relation 

^^^mA = S{t + tJ = R{t + tw,0)-R{t,tw) (19) 
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As far as the limit of R{t, tyj) when t ^ oo does not depend on the quantity defined in eq. 
(|19D goes to zero when t oo. In this way we can get the value of the integrated response by 
measuring the stress density as function of time. 
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Figure 2: On a logarithmic scale we see the response function 5 at F = 1.8 for 
different values of (128, 512, 2048, 8912 and 32768) as function of s = t/t^. 



130 at 



In our case (where we use the stress as a perturbation) the physical interpretation of the 
procedure is quite clear. We start by putting the systems in a box which is not cubic (because 
e 7^ 0), but two sides are slightly longer of the third. At time t^, we change the form of the 
box to a cubic one. In this way we deform the the system and we induce a stress which will 
be eventually decay. In the high temperature phase, where the system is liquid, the stress will 
disappear in a short time. On the contrary, in the glassy phase, we shall see that the stress 
remains for a much longer time (as expected IQ) and it shows an interesting aging behaviour. 

The choice of e is crucial. In principle its value should be infinitesimal. However the signal 
is proportional to e while the errors are e independent. The errors on the response function 
grow as e~^. On the other hand if we take a too large value of e we enter in a non linear region. 
We have taken data at e = .1 and e = .05 and we have seen that there are some non-linear 
effects. No non-linear effects have been detected at e = .02 and e = .01. All the data we 
present in this paper come from e = .01 and they are reasonable free of systematic effects. 

We plot the response function at F = 1.8 for A^ = 130 at different values of (128, 512, 
2048, 8912 and 31768) as function of s on a linear scale (|l]) and on a logarithmic scale @). 

As we can see the two regions FDT and aging are quite clear. It is also evident that S{go) 
is different from zero at this temperature (it is obviously zero in the high temperature phase, 
a residual dependence on t^. The short time region is shifted at smaller and smaller s when 
increases 

The crucial step would be now to plot the fluctuations C and the response S together. 
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The previous equation tell us that 

^-XiC) (20) 
SO that it is convenient to plot 5" versus C. The slope of the function S{C) is thus X{C). 




Figure 3: The response 5* as function of C at F = 1.8 for = 130 at different values of t^, 
(128, 512, 2048, 8912 and 32768). The two straight lines have slope 1 and .25. 

In the one step replica symmetry breaking scheme we expect that: 

S = mC for C<C(oo), 
S = C + {l-m)C{oo) for C < C(oo). (21) 

In fig. (^) we find the data for the response 5" as function of C at F = 1.8 for = 130 
at different values of ( 512, 2048, 8912 and 32768). We can see that apart from the data 
at tyj = 512 the data for higher are compatible, inside the errors and the predicted scaling 
works quite well. We can easily see a region where the FDT theorem is valid. Beyond this 
region the data are compatible with a linear dependance of S* on C as in the case of one 
step replica symmetry breaking. The straight line correspond to one step replica symmetry 
breaking with m = .25. This value of m is quite compatible with the value of m (m = .33 ±.04) 
found in a different simulation where a different quantity (i.e. the diffusion) was studied) 0. 

The numerical results are in very good agreement with the theoretical expectations and 
support the possibility of one step replica symmetry breaking for real glasses. 
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